Observational constraints and cosmological parameters 
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I discuss the extraction of cosmological parameter constraints from the recent WMAP 3-year 
data, both on its own and in combination with other data. The large degeneracies in the first 
year data can be largely broken with the third year data, giving much better parameter constraints 
from WMAP alone. The polarization constraint on the optical depth is crucial to obtain the main 
results, including ris < 1 in basic six-parameter models. Almost identical constraints can also be 
obtained using only temperature data with a prior on the optical depth. I discuss the modelling of 
secondaries when extracting parameter constraints, and show that the effect of CMB lensing is about 
as important eis SZ and slightly increases the inferred value of the spectral index. Constraints on 
correlated matter isocurvature modes are not radically better than before, and the data is consistent 
with a purely adiabatic spectrum. Combining WMAP 3-year data with data from the Lyman-a 
forest suggests somewhat higher values for as than from WMAP alone. 



I. INTRODUCTION 

The WMAP satellite has now provided three years 
worth of beautiful data giving a clear picture of the CMB 
temperature over most of the sky and allowing robust 
constraints on many cosmological parameters [3, 01- I 
shall briefly review parameter estimation methodology, 
then discuss in more detail topical issues relating to 
WMAP parameter constraints. I also discuss the utility 
of combining with other data, and show some example 
joint constraints obtained by combining with Lyman-a 
data. I shall assume familiarity with standard abbrevia- 
tions and cosmological parameters. 



II. PARAMETER ESTIMATION 

The CMB sky is expected (and observed) to be closely 
Gaussian, and noise properties can also often be ap- 
proximated as Gaussian. Given a set of parameters, 
linear-theory predictions for the CMB power spectrum 
can be computed quickly using numerical codes such as 
CAMB 3] and CMBFAST Q. Since the statistics of ob- 
served sky or power spectrum estimators are quite well 
understood, we can therefore compute the likelihood of 
a given set of parameters given the observed data. We 
then wish to calculate the posterior distribution of vari- 
ous parameters of interest, P(parameters|data), given a 
set of priors. In potentially large dimensional parame- 
ter spaces, the information in the posterior distribution 
is most easily extracted by using a set of samples from 
the distribution. For example an estimate of the poste- 
rior mean of a given parameter is given simply by the 
average of the value of that parameter in the set of sam- 
ples. More complicated quantities can also be computed 
from samples, for example the correlations of different 
parameters, marginalized two-dimensional distributions, 



etc. 

The methods almost universally adopted for sampling 
are variations of the Metropolis-Hastings Markov Chain 
Monte Carlo method. This method is explained in e.g. 
Refs. 0,001 3.nd consists of constructing a way of ex- 
ploring parameter space such that the value of the pa- 
rameters at any given step are a sample from the pos- 
terior distribution. There is considerable choice in the 
details of how these methods are implemented, and some 
discussion of suitable step-proposal distributions is given 
in Ref . |^ . A fairly complete implementation aimed 
at sampling and analysing non-pathological distributions 
(i.e. preferably unimodal and not too distorted) most 
common in cosmology is publicly available in the Cos- 
moMC^ package iflj used here; see also AnalyseThis p^ . 

Once a set of samples is available, almost every statisti- 
cal question of interest can be answered using estimators 
constructed from the samples. The probability density 
at a point in parameter space is given simply by the 
number density of samples. Marginalized probabilities 
in reduced dimensions are given simply by the number 
density of samples for the parameters in those reduced 
dimensions. Furthermore often new priors or data can 
be added very quickly using importance sampling: a way 
of re-weighting samples based on the difference of the 
new and old likelihoods — see Ref.i for further discus- 
sion. Chains for the WMAP 3-year analysis are available 
from the LAMBDA^ website, and are compatible with 
CosmoMC's GetDist program. 



III. WMAP 3- YEAR CONSTRAINTS 

It is well known from the WMAP first year data that 
with good measurements of only the first two acoustic 
peaks of the CMB temperature power spectrum there 
remains an important degeneracy between the spectral 
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FIG. 1: Constraints from WMAP 3-year temperature (points) 
and joint with polarization (68% and 95% contours) for a ba- 
sic six parameter ACDM model (no tensors). The points 
represent samples from the posterior distribution, and are 
coloured by the value of the optical depth r. Polarization 
constrains the optical depth, breaking the main flat-model 
degeneracy and suggesting < 1. 



index Ug, the optical depth r, the baryon density flhh^ 
and the amplitude of fluctuations parameterized by Ag or 
CTg. This degeneracy remains in the 3-year temperature 
data, as shown by the samples in Fig. 
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FIG. 2: Constraints from WMAP 3-year temperature (red), 
temperature and polarization (black), and temperature with 
a Gaussian prior on the optical depth r = 0.10 ± 0.03 (blue). 
The top six parameters have flat priors and are sampled using 
MCMC, the bottom six parameters are derived. 



A. Constraining the optical depth 

A large scale polarization signal on the CMB can be 
generated at reionization due to scattering of the CMB 
quadrupole. By measuring the large scale E-mode power 
spectrum, one can therefore constrain the optical depth. 
The 3-year WMAP analysis makes a heroic attempt to 
clean out dominating foregrounds to extract this sig- 
nal ^3 • Various consistency checks (for example B-mode 
power spectrum consistent with zero) suggest that this 
has been achieved to sufficient accuracy to reliably con- 
strain the optical depth. However the possibility of some 
unidentified foreground remains, and it is unclear how 
large the foreground-induced systematic error in the re- 
sult is (unlike for the temperature analysis, there is no 
marginalization over the foreground templates or fitting 
uncertainties). 

Using the WMAP 3-year result for the large scale po- 
larization the optical depth t can be constrained, break- 
ing the main flat-model degeneracy, and suggesting that 
rij < 1 in basic tensor-free models as shown in Fig. ^ 
In fact this is essentially the only information in the po- 
larization that effects basic cosmological model parame- 
ters. For example if we take a prior on the optical depth 
T — 0.10 ± 0.03 (consistent with the polarization result) 
essentially identical parameter constraints are obtained 



using only the temperature results as shown in Fig. |21 
The polarization spectrum on smaller scales may how- 
ever be useful for constraining extended models, for ex- 
ample with isocurvature modes. The improvement in 
parameter constraints from the better measurement of 
the temperature power spectrum with the 3-year data is 
relatively modest; using the same r prior with the first 
year temperature power spectrum gives an almost iden- 
tical constraint on Ug as using the 3-year data, though 
the matter densities are constrained rather better due to 
the better measurement of the second and third acoustic 
peaks. 

Different priors on the optical depth or reionization 
redshift can be applied very quickly to existing chains 
using importance sampling. For example, my personal 
prior is approximately a Gaussian with Zj-o = 8 ± 3 trun- 
cated to zero at z^c < 6, which is towards the lower end 
of the distribution obtained by WMAP (e.g. if there 
were an unidentified foreground remaining in the WMAP 
data, the value of the optical depth would be expected to 
come down further). In this case using WMAP temper- 
ature and polarization without tensors the < 1 result 
is very significant (n^ < 0.98 at 2-sigma) and the value 
of (Ts comes down a bit to 0.73 ± 0.05. 
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FIG. 3: The theoretical unlensed (black) and smoother lensed 
(red) CMB temperature power spectra (top) and the differ- 
ence (bottom blue) for a fiducial WMAP 3-year ACDM model 
with Us = 0.95, T = 0.09. Green lines are unlensed but in- 
clude the SZ spectrum from Ref. evaluated for fiducial 
parameters as used by the WMAP team (with marginalized 
amplitude) at 22.8GHz. At higher frequencies the SZ contri- 
bution is somewhat smaller. 




FIG. 4: Six parameter WMAP 3-yr constraints with no sec- 
ondary modelling (blue), marginalizing over SZ model am- 
plitude (red) and marginalizing over SZ and including CMB 
lensing (black). 



B. Effect of secondaries 

The WMAP 3-year parameter analysis accounts for an 
SZ contribution by marginalizing over the amplitude of 
an analytic model calculated for a fixed fiducial set of pa- 
rameters 2j . The effect of including SZ on the cosmolog- 
ical parameter constraints is not large, and comparable 
to that obtainable by changing priors in the base param- 
eters. Since the SZ adds power on small scales, adding a 
contribution means that Ug has to decrease slightly, shift- 
ing the posterior distribution further down from = 1. 
However when modelling secondaries, all relevant effects 
should be included. It is well known that CMB lensing 
also has a percent level effect on the relevant scales — 
see Fig. I^Cfor a recent detailed review of CMB lensing 
see Ref. ,1^). Lensing smooths out the acoustic peaks, 
and on WMAP scales the main effect is a reduction in 
the height of the third peak height, as well as smaller 
effects on the earlier peaks and troughs. The effect can 
be compensated by increasing rt^, so including lensing 
tends to increase Ug. The effect is shown in Fig 0] and 
is comparable magnitude to SZ, but has an opposite ef- 
fect on Tig. The best fit sample has a slightly better 
with lensing included, with Ax^ ^ 1. Since the effect 
of secondaries is small, they do not have a large effect 
on the conclusions of the WMAP 3-year analysis. How- 
ever if SZ is included it would seem sensible to also in- 
clude CMB lensing. The CMB lensing effect can be com- 
puted very accurately and is trivial to include, and the 
effect on WMAP parameters shown here is compatible 
with the forecast given in Ref. 14] . The fact that SZ 
and lensing have nearly opposite effects suggests that in 
practice rather accurate WMAP constraints on can 
be obtained by neglecting both. However the effect on 
ilbh^ is in the same direction, and including both shifts 
the mean up by nearly one sigma compared to including 
neither, giving Qbh'^ = 0.0228 ± 0.0008. 



C. Effect of likelihood estimation 

Though in principle the likelihood from a Gaussian 
CMB sky should be easy to calculate, in practice it can be 
rather difficult due to partial sky coverage, dataset size, 
non-uniform noise, beam uncertainty, etc. [J]. The 3-year 
WMAP likelihood approximation is much improved from 
the first-year one, including full covariance information 
on large scales where the C/ are most non-Gaussian and 
noise correlations are most important. Unless otherwise 
stated results quoted in this paper are computed using 
the likelihood code supplied by the WMAP team with 
default settings and no secondary modelling. However 
there are options in the likelihood code depending on how 
the effect of the beam uncertainty is modelled [13 • For 
basic six-parameter models, using the alternative options 
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FIG. 5: Constraints on the spectral index Us and as from 
WMAP only (green), SDSS Lyman-a (blue, Ref. [l3), and 
the joint constraint (red). Contours enclose 68% and 95% 
of the probability, and the model is ACDM with no tensors. 
Lyman alpha only contours are for fixed best-fit joint values 
of the other parameters and dependent upon this prior: they 
give a visual clue to the direction and amount by which the 
data pull the joint constraints, but the absolute position of 
the contours is fairly meaningless. 



of Ref. ^3 tends to shift n^. to slightly higher values'^, in- 
creasing Us by ^ 0.006. This suggests an estimate of the 
systematic error due to small scale likelihood modelling 
of 0(0.01) on Us, comparable to the effect of secondaries 
and choice of variables with flat priors. 



IV. COMBINING DATA 

WMAP alone provides good constraints on many pa- 
rameters, especially those that are directly related to the 
power spectrum. For example the ratio of the sound hori- 
zon at last scattering to the angular diameter distance 
to last scattering is constrained to ~ 1%. Other com- 
binations of parameter are also measured very well, for 
example WMAP alone gives the linear theory constraints 

V0.238j V0J3/ V 0-689 J ~ ' ^ ^ 

(Llfl) f M' 7ity°°" = i± 0.04.(2, 
V 0.689 J V0.73 J \0.238j ^ ' 
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FIG. 6; The 95%-confidence regions for the CMB polariza- 
tion power spectra given WMAP 3-yr, Boomerang and 
2dF 111 data with an HST prior h = 0.72 ±0.08. The param- 
eters varied are those for a constant spectral index ACDM 
model with insignificant neutrino mass and sharp reioniza- 
tion. The two tensor results are the B-mode power spectra 
expected from for a scale-invariant tensor-mode power spec- 
trum with At = 4r X 10~^ and two possible values of r. The 
effect of tensors on the cosmological parameter constraints 
was neglected, and a prior z^c > 6 was assumed. This is an 
updated version of the figure in Ref. 13]. 



However the third combination of these parameters is 
only constrained to 20%. Any data that can improve 
on these constraints in any direction can be used to im- 
prove parameter constraints. In addition external data 
provides an important consistency check, and any devia- 
tions from consistency may indicate departures from the 
assumed cosmological model. 

The WMAP results are now so good that for simple 
models the small scale C; can be predicted rather ac- 
curately from the WMAP data on larger scales. Data 
from other CMB experiments currently give relatively 
large errors on the smaller scale acoustic peaks, and 
the improvement to parameter constraints from includ- 
ing other CMB data on linear scales is now rather mod- 
est. To significantly improve parameter constraints from 
3-year WMAP a tight constraint on the third acoustic 
peak would be very useful. As shown in Fig. |S1 small 
scale polarization observations have to be very sensitive 
to compete with the prediction that can be made from 
current data assuming a basic cosmological model, and 
it will be a while before small scale CMB i5-polarization 
helps at all with constraining parameters in these models. 
However consistency with the predictions is an important 
consistency check. 
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FIG. 7: Allowing tensors, samples from WMAP 3-year tem- 
perature and polarization (points), compared to the joint con- 
straint with SDSS Lyman-a data from Ref. ^ (68% and 95% 
contours). 



The WMAP team give constraints with various combi- 
nations of data as discussed in Ref. 2] . Of particular note 
is the significant increase in the expected value of erg when 
including data from weak lensing. Here as an example 
I consider combination with Lyman-a data from SDSS, 
Refs. 17] . This serves to measure the matter power spec- 
trum scales much smaller than those directly probed by 
the CMB and hence has the potential to beak degenera- 
cies involving and erg- The data from Ref. fliil are less 
constraining that the SDSS results, so I concentrate on 
the latter here. 

Figure |S1 shows the joint constrain from WMAP and 
Lyman-a without tensors. As with weak lensing, adding 
Lyman alpha increases the values of (Tg from the value 
preferred by WMAP alone, with the SDSS joint con- 
straint giving (Tg = 0.86 ± 0.03. The joint constraint is 
actually otherwise fairly consistent with WMAP alone, 
for example the best fit sample has r — 0.09, rts — 0.97, 
though it does favour higher values of fim ~ 0.3. It is 
consistent with Us < I at two sigma with no running of 
the spectral index. 



V. EXTENDED MODELS 

One of the most interesting extensions of the basic 6- 
parameter ACDM model is one allowing for a contribu- 
tion from tensors modes as predicted by some models of 
inflation. FigureHshows constraints on the tensor-scalar 
ratio r (defined in terms of the initial power spectra as in 
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FIG. 8: Constraints on the correlated matter isocurvature 
mode ratio B as defined in Ref. f2 1| using WMAP 3-year data 
combined with small scale CMB |23,|23 and 2dF data [l^. 



the WMAP papers) from WMAP alone'* and when com- 
bined with Lyman-a. Clearly there is no evidence for a 
tensor contribution, however allowing for tensors does al- 
low the possibility of > 1: if tensors add power to the 
low-^ multipoles, the spectral index can be bluer and still 
maintain the correct overall ratio of small and large scale 
power. Many future observations will be aimed at detect- 
ing the i?-mode polarization signal from tensor modes, 
and Fig. |H| shows the expected spectrum for a couple of 
possible amplitudes in comparison to the predicted sig- 
nal from lensed scalar modes. Further discussions of in- 
flationaryconstraints from WMAP 3-year data are given 
in Refs. Il[ll|2|. 

Another interesting possibility is an isocurvature mode 
contribution. Probably the simplest possibility is a to- 
tally correlated contribution from matter isocurvature 
modes, for example from a curvaton model of inflation. 
Here the spectral index is fixed to that of the adiabatic 
modes, and there is only one extra parameter, the ratio 
of isocurvature modes (which can be negative for anti- 
correlated modes). The updated constraint from Ref. [^J 
is shown in Fig.|Sl corresponding to a 95% confidence con- 
straint —0.42 < B < 0.25. As with older data there is 
no evidence at all for isocurvature modes, which should 
not be surprising given that a simple adiabatic model 
fits the data well. The WMAP 3-year data does not rad- 
ically improve the constraint because the main effect is on 
large scales where the first year data was already cosmic 
variance limited in the temperature. However degener- 
acy breaking does help, as does improved sensitivity on 
the polarization. I expect that more general isocurvature 
analyses (e.g. updating Ref. 24]) should have a similar 



* Note that the contours in the corresponding Fig. 14 of Ref. 
were wrong in version 1 of that paper. 
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conclusion: no evidence for isocurvature modes, though 
a significant fraction stiU aUowed. 

VI. CONCLUSIONS 

The WMAP 3-year results appear to show a remark- 
able agreement with a simple six parameter ACDM 
model. Numerous extended models have been discussed 
by the WMAP team P|, though none fit the data signif- 
icantly better. There is some tension by between the as 
values favoured by WMAP and weak lensing and Lyman- 
alpha data, and this deserves further investigation. The 
WMAP data alone suggest that either Ug < 1 or there 
is something else interesting, for example a significant 
tensor mode component. 

Future CMB observations that can provide an accurate 
measurement of the optical depth or third acoustic peak 
are needed to constrain parameters significantly better 
with CMB alone. Precision measurements of the temper- 
ature spectrum to small scales should allow the spectral 
index to be determined independently of the polariza- 
tion constraint on the optical depth. Current and future 
small scale polarization observations are unlikely to im- 
prove basic parameter constraints significantly, though 
they can be a good test of alternative models. Detecting 
large scale S-polarization from gravitational waves is of 
course the main goal of many future missions, and would 

^ jhttp : //lambda .gsfc .na sa.gov/ 1 



be a powerful way to constrain models of inflation. 

CMB observations are in principle clean way to con- 
strain early universe perturbations and many cosmolog- 
ical parameters. However they are rather insensitive to 
other parameters, for example the late time evolution of 
the dark energy. To constrain these parameters com- 
plementary data are needed. Consistency with different 
data sets is a good check on the assumed cosmological 
model. If consistent other data sets can also be used to 
break partial degeneracies in the CMB constraints, giv- 
ing better joint parameter estimates. Many additional 
data sets not considered here are analysed in combina- 
tion with WMAP 3-year data in Ref.^. 

Acknowledgements 

I thank the organisers for inviting me to give this talk 
and an excellent conference. I thank Matteo Viel, George 
Efstathiou, Hiranya Peiris and fellow conferees for use- 
ful discussion. I thank K. Abazajian and J. Lesgour- 
gues for making their Lyman-a CosmoMC modules pub- 
lic. I acknowledge the use of the Legacy Archive for Mi- 
crowave Background Data Analysis^ (LAMBDA). Sup- 
port for LAMBDA is provided by the NASA Ofhce of 
Space Science. I thank CITA for use of their Beowulf 
computer. I am supported by a PPARC Advanced Fel- 
lowship. 



[1] G. Hinshaw et al. (2006), astro-ph/0603451" 

[2] D. N. Spergel et al. (2006), astro-ph/0603449 

[3] A. Lewis, A. Challinor, and A. Lasenby, Astrophys. J. 

538, 473 (2000), astro-ph/9911177 
[4] U. Seljak and M. Zaidarriaga, Astrophys. J. 469, 437 

(1996), astro- ph/9603033^ 

[5] R. M. Neal (1993), http : TTc osmologist . inf o/Neal93] 
[6] N. Christensen, R. Meyer, L. Knox, and B. Luey, Class. 

Quant. Grav. 18, 2677 (2001), astro-ph/0103134 
[7] A. Lewis and S. Bridle, Phys. Rev. D66, 103511 (2002), 

|astro-ph/0205436| 
[8] A. Lewis Mid Bridle, 

"http : //cosmologist . inf o/notes/cosmomc .ps .gz" 
[9] J. Dunkley, M. Bucher, P. G. Ferreira, K. Moodley, and 

C. Skordis, Mon. Not. Roy. Astron. Soc. 356, 925 (2005), 

astro-ph/0405462 
[10] M. Doran and C. M. Mueller, JCAP 0409, 003 (2004), 

astro-ph/0311311 
[11] L. Page et al. (2006), p[st ro-ph/0603450] 
[12] E. Komatsu and U. Seljak, Mon. Not. Roy. Astron. Soc. 

336, 1256 (2002), astro-ph/0205468 
[13] A. Lewis and A. Challinor (2006), astro-ph/0601594 



[14] A. Lewis, Phys. Rev. D71, 083008 (2005), 

astro-ph/0502469 
[15] H. Peiris and R. Easther (2006), astro-ph/0603587 
[16] M. Viel, M. G. Haehnelt, and V. Springel, Mon. Not. 

Roy. Astron. Soc. 354, 684 (2004), astro-ph/0404600 
[17] P. McDonald et al., Astrophys. J. 635, 761 (2005), 

astro-ph/0407377 

[18] W. C. Jones et al. (2005), astro- ph/0507494 

[19] W. J. Percival et al., MNRAS 327, 1297 (2001), 

astro-ph/0105252 ^ 

[20] L. Alabidi and D. H. Lyth (2006), astro-ph/0603539" 
[21] C. Gordon and A. Lewis, Phys. Rev. D67, 123513 (2003), 

astro-ph/0212248 
[22] A. C. S. Readhead et al., Astrophys. J. 609, 498 (2004), 

astro-ph/0402359 
[23] C.-l. Kuo et al. (ACBAR), Astrophys. J. 600, 32 (2004), 

astro-ph/0212289 
[24] M. Bucher, J. Dunkley, P. G. Ferreira, K. Moodley, 

and C. Skordis, Phys. Rev. Lett. 93, 081301 (2004), 

astro-ph/0401417j 



